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Electron-hole generations: A numerical approach to interacting fermion systems 
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A new approach, motivated by Fock space localization, for constructing a reduced many-particle 
Hilbert space is proposed and tested. The self-consistent Hartree-Fock (SCHF) approach is used to 
obtain a single-electron basis from which the many-particle Hilbert space is constructed. For a given 
size of the truncated many particle Hilbert space only states with the lowest number of particle-hole 
excitations are retained and exactly diagonalized. This method is shown to be more accurate than 
previous truncation methods, while there is no additional computational complexity. 
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There has been a growing interest in the physics of 
low temperature disordered many-particle systems such 
as the two dimensional "metal- insulator" transition 
electron dephasing due to interactions Q and transport 
through quantum dots || . Unfortunately these problems 
are very difficult to treat analytically since neither disor- 
der nor interactions may be considered as a small pertur- 
bation. An exact numerical treatment is limited to small 
systems due to the large size of the many-particle Hilbert 
space. For a system of m single particle orbitals and n 
electrons, the many-particle Hilbert space size M = (™). 

One popular approximation method for low-lying 
eigenstates (especially the ground state) is to reduce 
the size of the many-particle Hilbert space, i.e., to limit 
the number of many-particle basis states spanning the 
Hilbert space. One strategy is to include in the many- 
particle basis states only single-electron states which 
are lower in their non-interacting energy than a certain 
threshold. This will limit the number of single-electron 
states ttir < m, reducing the number of many-particle 
states to Mr = (™ R ). This procedure works well in the 
strongly localized regime H||, but fails for weaker dis- 
order. As demonstrated bellow, a similar strategy using 
Hartree-Fock single-electron states, works better in the 
diffusive regime. 

One might also reduce M by imposing a global con- 
straint on the many-particle states. For example, the 
many-particle basis states may be composed of Mr < M 
Fock determinants of single-electron states for which the 
total non-interacting energy or the expectation value of 
the full interacting Hamiltonian is smaller than some 
threshold |J . An even better procedure is to use Hartree- 
Fock single-electron states to build the many-particle ba- 
sis, while picking the Mr states with the lowest values 
of the expectation value of the interacting Hamiltonian 
|?||| . All these methods suffer from the drawback of forc- 
ing one to go over all M basis states in order to chose Mr 
of them according to the global criteria. Since M grows 
exponentially, this is not a practical solution. Neverthe- 
less, it is possible to overcome this snarl by sampling only 
part of all the M states, for example by a Monte Carlo 



procedure M. 

In this paper we propose a different way to trun- 
cate the many-particle basis which gives a more accurate 
evaluation of the low-lying eigenvalues and eigenvectors 
for the same truncated basis size Mr than the previ- 
ous methods. Moreover, the truncation criteria is not 
global and does not require evaluating some criteria for 
all the many-particle basis states, i.e., the truncation has 
the same computational complexity as the single-electron 
criteria truncation but nevertheless produces more accu- 
rate results than the complicated global criteria meth- 
ods. Our truncation method is based on the fact that 
the interaction operator is a two-body operator. Thus, 
only many-particle states which differ by at most two 
e-h (electron-hole) pairs can be coupled by the interac- 
tion operator. This fact is well known in the context 
of nuclear physics |^,|o|, and has recently been brought 
to the forefront of disordered interacting systems by the 
concept of Fock space localization by Altshuler, Gefen, 
Kamanev and Levitov [jll| . We therefore propose to con- 
struct the many-particle basis in the following way: First 
a many-particle basis state is constructed as a determi- 
nant from the n lowest-lying Hartree-Fock single-electron 
states. Then all the many-particle basis states differing 
from the first state by up to k e-h pairs are included, 
resulting in a truncated basis of size 



Mi 
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From the point of view of the Fock space localiza- 
tion postulate [|ll| this truncation method has the po- 
tential for being accurate for finite systems provided k 
is larger the Fock space localization length £p, since 
the average contribution of a state in the fc-th gener- 
ation to the exact eigenstate of the system is propor- 
tional to exp(— fc/£p). Thus, the contribution of the 
state of the A:-th e-h generation to the exact eigenstate 
is proportional to (™ _n ) (£) exp(— fc/£ir), i.e., as long as 
1/£f > ln((m— n)n) the weight of each consecutive gener- 
ation falls of exponentially. Although there is an ongoing 
debate regarding the critical energy at which the transi- 



tion occurs and whether the transition is abrupt or more 
gradual ]l^-pj| , it is nevertheless clear that for low- lying 
excitations the contribution of large k states will fall off 
strongly. 

We illustrate the method for an interacting disordered 
many-particle two dimensional tight-binding Hamilto- 
nian: 

H = H o + H mi 

H = Y^ e k,j n k,j - V ^2[ a k,j+l a k,j + a k+l,j a k,j + h - c -} 
k,j k,j 

Hmt = U 2J n ktj n k ±ij±i, (2) 

where {k, j} denotes a lattice site, a\. ^ is an electron cre- 
ation operator, the number operator is rik,j — a \j a k,j, 
£k.j is the site energy chosen randomly between — W/2 
and W/2 with uniform probability and V is a constant 
hopping matrix element. Only nearest neighbor interac- 
tions U are taken into account. This is reasonable in the 
metallic regime where screening is good. 

Systems composed of a n — 4, 6, 10 electrons residing 
on m — 16, 15, 24 sites of 4 x 4, 5 x 3, 6 x 4 lattice (i.e., 
M = 1820,5005,1961256), were considered. The latter 
size is at the limit of the possibility for exact diagonal- 
ization using the Lanczos method for current computer 
technology while the former sizes are typical for standard 
exact diagonalization procedures. Disorder strength of 
W = 8V for the 4 x 4 and 5 x 3 lattices and W = 5V for 
the 6x4 lattice were chosen so that the single particle 
localization length is larger than the system size, while 
the single-electron states follow the Gaussian orthogonal 
ensemble (GOE). 

We carry out exact diagonalization of the many- 
particle Hamiltonian using standard procedures for the 
smaller systems and the Lanczos method for the larger 
ones, obtaining the exact many-particle eigenvalues E a 
and eigenvectors \a), where a = 0,1,2 for the ground 
state, the first and second excitation. 

In order to obtain the single-electron states we diag- 
onalize the single electron part Hq in the Hamiltonian 
given in Eq. (^|), obtaining the single-electron eigen- 
vectors \ipi) and eigenvalues E{. For the Hartree-Fock 
single-electron wavefunctions we solve the self consistent 
Hartree-Fock (SCHF) Hamiltonian: 

H HF = H + ^ n k,jU(rik±i,j±i}o (3) 

~ ^2 4j a fc±ij±i^( a l±ij±i a fcj)o, 

k,j 

where (. . .)o denotes an average on the ground state 
which is calculated self-consistently. The SCHF single- 
electron eigenvectors \tpi) and eigenvalues £i are obtained 
through a self consistent diagonalization of the Hartree- 
Fock Hamiltonian. 



Any many-particle wavefunction \a) may be expressed 
in terms of a combination of Fock determinants composed 
of the single-electron (or Hartree-Fock single-electron) 
eigenvectors in the following way: 

i«> = c iui2,..,i n 4 n ■ ■ ■ 4Ai i vac >' ( 4 ) 

where the sum is over all the possible combinations of N 
states out of the M possible ones for which m > i n > 
i n -i > . . . > %2 > it > 1, cj is the creation operator of 
the i-th. single-electron (or Hartree-Fock single-electron) 
state, vac) is the vacuum state. 

Let us now define more precisely the various trunca- 
tion methods which we have discussed above. The first 
method (which we shall name single-electron energy trun- 
cation) is to include in the many-particle basis states only 
single-electron states for which £j < er. This method will 
reduce the number of single-electron states to ttir < m. 
Usually £r is chosen according to the maximum m# de- 
sired. The Mr = (™ R ) many-particle basis will include 
the states: 

|/3) = c] n ... C J 2C l i |vac) ) (5) 

for which ttlr > i n > i n —i > . . . > ii > i\ > 1. The vari- 
ant of this method using the single-electron Hartree-Fock 
states for which C,i < Cr will be called the single-electron 
Hartree-Fock energy truncation. The basis in this case 
composed of the states described in Eq. (||) in which 
now c| is the creation operator or the i-th Hartree-Fock 
single-electron state. 

The truncated Hamiltonian then may be diagonalized 
using standard diagonalization procedures as long as Mr 
is not too big (Mr ~ 10000). For larger sizes one can 
use the Lanczos method as long as the matrix remains 
sparse enough. The sparseness of the many-particle ma- 
trix is roughly S ~ n 2 (mR — n) 2 / (™ fi ). Thus, as long 
as m,R is not to close to n the Lanczos method may be 
useful in obtaining the low-lying excitations. 

For the methods based on truncation according to a 
global criteria, one (in principal) should construct all 
possible basis states |/3) for which m > i n > i n -i > 
. . . > 12 > ii > 1- Then only states which correspond to 
some global criteria, for example (j3\H\(3) < Er, are to 
be included in the basis. Again this method may use ei- 
ther \ifi) or as the single-electron states from which 
the many-particle basis states |/3) are composed. This 
method will be named global energy truncation. In this 
case the sparseness of the matrix depends on the global 
criteria, and for a small Mr the matrix is not necessarily 
sparse JT). 

In the e-h generation truncation method the basis 
states are composed as follows: A basis state |/3 ) = 
c\ L . . . c\c\ |vac) is chosen. Then all states which may be 



created from this state by up to k applications of an e-h 
pair creation operator are added to the basis, i.e., 

= n* =1 Ifl,), (6) 

where m > l± ^ I2 ^ . ■ ■ h > n and n > i\ ^ 
i2 7^ . . . ik > 1- The total number of basis states 
for up to k e-h generations is given in Eq. (|l|). The 
sparseness of a specific e-h generation j is given by 
Sj ~ (n - j) 2 (m -n- j) 2 / (™T 2 ") (?_ 2 ), and the many- 
particle matrix sparseness will be determined mainly by 
the higher generations. Thus, the higher the number of 
e-h generations included, the sparser the matrix becomes. 

In the following we shall compare the accuracy of the 
different methods. First we compare the accuracy of the 
single-electron energy truncation to the e-h generation 
truncation method by applying both of them to the 5x3 
system. In Fig. |l| a comparison of the overlap between 
the approximated wave function \a') obtained using the 
truncated basis {|/3}} and the exact wave function \a) 
is presented. The results were averaged over 10 differ- 
ent realizations of disorder. For the single-electron en- 
ergy truncation we have chosen er for each realization 
such that values of rriR = 6, 7, 8, 9, 10, 11, 12, 13, 14 and 
Mr = 1,7,28,84,210,462,924,1716,3003 correspond- 
ingly are obtained. This is compare with the results 
of the generation truncation of k = 0,1,2,3,4,5 and 
Mr = 1, 55, 595, 2275, 4165. Both cases of single-electron 
energy truncation (Fig. |l|a) and single-electron Hartree- 
Fock energy truncation (Fig. |l|b) were considered. It 
can be seen that for roughly the same size of the trun- 
cated many-particle basis Mr the generation truncation 
is better by up to an order of magnitude for either the 
single-electron or the Hartree-Fock single-electron states. 
It is also interesting to note that using Hartree-Fock 
states does not always improve the accuracy. While for 
the ground state the Hartree-Fock single-electron states 
give a better result in any truncation method than the 
non-interacting single-electron states, for the first excited 
state this is so only as long as Mr is small compared to 
M. 

In Fig. H we check the accuracy of the differ- 
ent truncation methods for the eigenvalues of a much 
larger 6x4 system. The results were averaged over 
100 realizations. For the single-electron energy trun- 
cation we have chosen sr for each realization such 
that values of m R = 10, 11, 12, 13, 14, 15, 16, 17, 18 and 
Mr = 1, 11, 66, 286, 1001, 3003, 8008, 19448, 43758 corre- 
spondingly are obtained. This is compared with the re- 
sults of the e-h generation truncation of k = 0, 1, 2, 3 and 
Mr = 1,141,4236,47916. The most obvious observa- 
tion is that using Hartree-Fock states in any truncation 
scheme is much better than using non-interacting single- 
electron states. This stems from the fact that for this 
large system Mr <C M. Again, for approximately the 
same many-particle basis size Mr the generation trunca- 



tion method is more accurate. 

A more stringent test for the generation truncation 
method is its comparison to the global energy truncation. 
In Fig. H the accuracy of the three lowest eigenvalues for 
the generation truncation and global energy truncation 
method are compared for a 6 x 4 system. The energy 
generation truncation after k = 1,2,3 generations (i.e., 
Mr = 141,4236,47916) is compared to a global energy 
truncation of the same basis size when the single level 
states in both cases are composed of Hartree-Fock states. 
For both methods the accuracy is better for the ground 
state than for higher excitations. In all cases the gen- 
eration truncation performs better at a given Mr than 
the global energy truncation. Moreover, the advantage 
of the generation truncation is larger for the lower values 

Of Mr. 

Let us now examine whether one can formulate a quan- 
titative estimation for the contribution of different e-h 
generations. The weight of a given Fock determinant 
in the exact eigenvector is \Ci lt i 2w .. t i n | 2 (see Eq. (||)). 
By averaging over all Fock determinants belonging to a 
certain e-h generation k one obtains the average weight 
of a state in the k-th generation (C?m) depicted in Fig. 
f§ The average weight roughly falls off exponentially 
as a function of k, (Cf k \) oc exp(— k/^p), where £f is 
the Fock space localization length. We find that £p de- 
pends on the interaction strength £,f(U = 2V) ~ 0.4 and 
£f{U = W) ~ 0.7 and has also a weaker dependence 
on m and n. This exponential behavior is in line with 
what is expected by the Fock space localization scenario 
. The number of states in the fc-th e-h generation is 
equal to (™~") (J?), which for k <C n, m is proportional to 
Um — n)n) k exp(— (m/2(m — n)n)k 2 ) / (k\) 2 (see inset Fig. 
Q). Thus, as long as 1/^f > hi[(m - n)n/(k + l) 2 ] the 
weight of each consecutive generation beyond the k-th 
generation becomes exponentially smaller. 

In conclusion, using the ideas of Fock space localization 
we have proposed and tested a new truncation method for 
reducing the size of the many-particle Hilbert space. The 
method is based on retaining states with the lowest num- 
ber of electron hole excitations above the SCHF ground 
states. This method is shown to be more accurate than 
other widely used truncation methods, while requiring 
no additional computational complexity. Refinements of 
the e-h generation truncation methods (e.g., adding a 
single-electron energy truncation condition) merit further 
study. 

I would like to thank Y. Gefen and B. L. Altshuler for 
useful discussions and the Israel Science foundation for 
financial support. 
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FIG. 1. The discrepancy between the exact eigenstate \a) and the one \ai) obtained via a truncated many-particle basis for 
n = 6 electron on m = 15 sites of a 5 x 3 lattice (i.e., M = 5005) for U = 2V. The discrepancy is presented as a function of 
the size of the truncated basis Mr, the state ( the ground state |0) (circles); first excited state |1) (squares)) and the method 
of truncation, (a) The single-electron energy truncation (dashed line), generation truncation method (full line), (b) As in (a) 
but the single-electron states are replaced by Hartree-Fock states. 




FIG. 2. The discrepancy between the exact eigenvalue E a and the one E a l obtained via a truncated many-particle basis for 
n — 10 electron on m = 24 sites of a 6 x 4 lattice (i.e., M = 1961256) for U = 3V. The discrepancy is presented as a function 
of the size of the truncated basis Mr, the energies ( the ground state Eo (circles); first excited state E\ (squares)) and the 
method of truncation, (a) The single-electron energy truncation (dashed line), generation truncation method (full line), (b) 
As in (a) but the single-electron states are replaced by Hartree-Fock states. 
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FIG. 3. The discrepancy between the exact eigenvalue E a and the one E a l obtained via a truncated many-particle basis for 
n = 10 electron on m = 24 sites of a 6 x 4 lattice (i.e., M = 1961256) for U = 3V. The discrepancy is presented as a function of 
the size of the truncated basis Mr, the energies ( the ground state Eq (circles); first excited state E\ (squares); second excited 
state E2 (diamonds)) and the method of truncation. The dashed line corresponds to the global energy truncation, while the 
full line corresponds to the generation truncation. In both cases Hartree-Fock single-electron states were used. 
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FIG. 4. The average weight of the exact ground state jO) (circles) and first excited state |1) (squares) per a k-th generation 
(C? fc )) for n = 6 electron on a 5 x 3 (m = 15) lattice at U = 2V (white symbols) and U = 4V (filled symbols), and for n = 4 
electron on a 4 x 4 (m = 16) lattice at U = 2V (grey symbols). All cases were averaged over 10 realization of disorder. A best 
fit to an exponential exp( — fc/^r) is indicated by the straight lines. Inset: the number of states per a k-th generation (™~") Q) 
is indicated by the diamonds (white symbols n = 6, m = 15, grey symbols n — 4, m = 16). The dashed line corresponds to the 
approximation ((m — n)m) k exp( — (m/2(m — n)n)k 2 )/(k\) 2 . 



